Marginal Fermi liquid resonance induced by a quantum magnetic impurity in d-wave 

superconductors 
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We consider a model of an Anderson impurity embedded in a d x 2_ y 2-wa,ve superconducting state 
to describe the low-energy excitations of cuprate superconductors doped with a small amount of 
magnetic impurities. Due to the Dirac-like energy dispersion, a sharp localized resonance above 
the Fermi energy, showing a marginal Fermi liquid behavior (lu\iilu as ui — > 0), is predicted for the 
impurity states. The same logarithmic dependence of self-energy and a linear frequency dependence 
of the relaxation rate are also derived for the conduction electrons, characterizing a new univer- 
sality class for the strong coupling fixed point. At the resonant energies, the spatial distribution 
of the electron density of states around the magnetic impurity is calculated, to be confronted with 
measurements of the scanning tunneling microscopy on Bi2Sr2Ca(Cui_ a; Ni a; )Og+5. 



PACS numbers:75.20.Hr,71.10.Hf,71.27.+a,71.55, 

It has been for a long time that non-magnetic and 
magnetic impurities are exploited to elucidate the mi- 
croscopic nature of the superconducting state. In the 
high-T c cuprates, substitution by divalent metals (Zn 
and Ni) for Cu in the Cu02 plane offers a particularly im- 
portant way of introducing such impurities, as they pre- 
serve the doping level and introduce only minimal struc- 
tural disorder. Recently, scanning tunneling microscopy 
(STM) has been further developed to probe the quasi- 
particle scattering states around a single Zn impurity in 
Bi2Sr2Ca(Cui_ x Zn a: )08+5 with a high spatial and en- 
ergy resolution by Pan et al. Q. In the obtained STM 
spectra, an intense zero-bias quasiparticle scattering reso- 
nance is found at the Zn sites, and the spatial dependence 
of the density of states (DOS) in the vicinity of the im- 
purities reveals a fourfold symmetry which characterizes 
d x 2_ y 2 wave superconductivity (dSC) In fact, the 
existence of a non-magnetic impurity induced resonant 
state in dSC was predicted theoretically by Balatsky et 
al. |-|| earlier. 

Here, we consider the scattering effects of magnetic im- 
purities on dSC. In Ni-doped cuprates, in order to main- 
tain the valence of Cu 2+ ions (3d 9 , S=l/2), the Ni 2+ 
ions have a configuration (3d 8 , S=l), and strong antifer- 
romagnetic exchange couplings with the neighboring Cu 
sites lead to a residual S=l/2 on the Ni site, acting as 
a localized magnetic spin weakly coupled to its environ- 
ment through exchange interactions [||. Although there 
are some theoretical studies on the quasiparticle states 
around such a magnetic impurity in dSC state 
most of them treated a classical magnetic impurity. It 
is thus interesting and timing to study how a quantum 
magnetic impurity affects the quasiparticle states in opti- 
mally doped Bi2Sr2Ca(Cui_ 2: Ni 2 ;)08+5 superconductors. 
Since the quantum fluctuations of the internal degrees of 
freedom of the magnetic impurity play an important role 
in the ordinary quantum impurity problems, one can thus 



expect that the effects of a quantum magnetic impurity 
are significantly different from those of a classical mag- 
netic impurity. Within the slave boson mean field (MF) 
theory, we predict a sharp resonance above the Fermi level 
and a marginal Fermi liquid (MFL) behavior for both im- 
purity and surrounding conduction electrons. Moreover, 
we explicitly calculate the spatial distribution of conduc- 
tion electron DOS to be compared with STM measure- 
ments. 

In this Letter, we assume that a BCS-type weak cou- 
pling theory is applicable as a phenomenological model 
for High-T c optimally doped superconductors though the 
underlying mechanisms are different. We also assume the 
magnetic impurities are described by the Anderson model 
with a strong Hubbard repulsion. When the correlations 
between the magnetic impurities on different sites are ig- 
nored, the model Hamiltonian is given by 
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where £k = h 2 k 2 / (2m) — is the dispersion of the con- 
duction electrons, Ak = Ao cos 2ip is the dSC order pa- 
rameter, and Ao is the gap amplitude. When Nambu 
spinors are introduced 
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we can simplify the model Hamiltonian in a matrix form 

U = £ ^k ( £ k^ + A k a x ) ^ k + V Y ($k f*£ + h.c^J 



+ (td + ^((p+a^+l) - ^($ + V - if, 
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where <j z and a x are Pauli matrices. Using the method of 
equations of motion, the generalized T-matrix is derived 
f (to + i0+) = Va z G d (w + i0+) Va z , where G£ (iu n ) = 
(ioj n — e^(T z — A^(T X ) 1 is the unperturbed Green func- 
tion (GF) of the conduction electrons. At zero tempera- 
ture, analytical continuation is used to calculate the per- 
turbed GF through the GF of the impurity: G (r, r'; u>) = 
G° (r - r', u) + G° (r, u) f (w) G° (-r', w). The local 
DOS of the conduction electrons around the magnetic 
impurity is thus given by N (r,u>) = — — ImGn (r, r;u>), 
and the relaxation rate for the conduction electrons is 
also deduced from ImTn(o; + i0 + ). 

When we take the infinite U limit, the impurity op- 
erator is expressed as (p + = ( f?b, f[b + ) in the slave- 
boson representation where the fermion f a and 
the boson b describe the singly occupied and hole states, 
respectively. The constraint b + b + fa fa = 1 nas to 

a 

be imposed. When a MF approximation is applied, the 
boson operators b and b + are replaced by a c-number 60, 
and the constraint is satisfied by introducing a chemical 
potential Ao- Therefore, the MF Hamiltonian is written 
as 



H ™f = ( £ k^ + Ak^) ^k + V (i>£ 
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where (f> + = ( f^~, f± ) denotes the Nambu spinors of the 
impurity quasiparticles and the renormalized parameters 
(d = f-d + A and V = b V. 

Using the standard techniques we find Gt(iu) n ) — 
1 _1 

j n — e d o~ z — £/ (iu> n ) , where the self-energy of the 
impurity becomes diagonal, 



Solving these equations yields &o and Ao for given param- 
eters W, A , T = ttN f V 2 , and e d , where N F is the DOS 
at the Fermi surface. In the following, we will choose Ao 
as the energy unit, W/A = 20, and r/A = 0.2. 

In the present model the local DOS of quasiparticles 
near the Fermi surface goes to zero linearly, so the usual 
logarithmic Kondo singularity in the scattering matrix 
of the magnetic moment with conduction electrons is 
thus absent. When e d is less than a threshold value, 6§ 
is zero, leading to a decoupled free local magnetic mo- 
ment, namely, no Kondo effect occurs. However, above 
the threshold value of e d , b^ rises steeply, and then sat- 
urates quickly. The usual broad mixed valence regime 
shrinks to a very narrow regime. In Fig. 1, the ground 
state phase diagram is calculated in the e d — V plane. 
For a given value of e d , there will be a phase transition 
from the decoupled free spin to the mixed valence, and 
finally a crossover to the strong coupling regime. The fi- 
nite threshold value of the phase transition is delineated 
by the solid line between area I and II, and turns out to 
be linear in je^l /A , approximately. 

In the mixed valence and strong coupling regimes, the 
impurity DOS versus e d is plotted in Fig. 2. A sharp 
local resonance always appears above the Fermi energy 
for each value of e d , while the corresponding DOS for 
lu < is broad and small. This is one of the most impor- 
tant differences between the magnetic and non-magnetic 
impurities scatterings, as the localized resonance always 
occurs below the Fermi energy for the repulsive potential 
scattering in the latter case . To the logarithmic accu- 
racy, the zero of the denominator of Gf{iuj n ) is given by 
Q = Q'-iQ", and 
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because of the inversion symmetry in the k summation. 

At T = 0, the ground-state energy change due to im- 
purity is: 



Simp = £d+A (&o - l)-- / dculn 



lu 2 (1 + a (lu)) 2 + el 
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where W is the band width and a (u>) = Z^+e' 2 a 
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The saddle-point equations are derived as 
w 



1 f 2uj 2 {l + a(Lo))a(u)) 

TT 





w 



bl 



cu 2 (l + a(Lu)) 2 + e 2 
2i d 



lu 2 (1 + a (u)T + ? d 



(5) 
(6) 



where VL' represents the position of the quasiparticle res- 
onance, while Q" corresponds to its width or the inverse 
lifetime. If the self-energy Ej (u>) is expanded near the 
resonant energy, the impurity DOS can be approximately 

written in a Lorentzian form Ni m „ (lu) « --, — , n ,,^> ■ 

imp v I 7t (a;— f2') +(0") 

At the resonant energy u = ft' the height of the res- 



onance is 
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inversely proportional to the resonant 
energy. As e d — > 0, this resonance becomes arbitrarily 
sharp and close to the Fermi surface, but the DOS at 
the Fermi energy is always suppressed to zero for all val- 
ues of e d because of the imaginary part of the impurity 
self-energy. 

Actually, an analytic expression for the retarded self- 
energy of the magnetic impurity can be derived, 
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for e = |w|/Ao < 1. Here, K is the complete elliptic 
integral of the first kind. As u> — > 0, we have 



ReE/(u;) 



-68 



2r 



lu In • 



4A n 



(7) 



2 



i m ^/H--^(M + ^), (8) 

i.e. precisely the MFL behavior proposed by Varma et 
al. to describe the anomalous normal state properties 
of optimally doped cuprates O]. Within the T- matrix 
approximation, the self-energy of the conduction electron 
has exactly the same type of singular behavior. To our 
knowledge, it is the first time to obtain such a result. 

Earlier, the Kondo effect in "gapless" fermion systems 
(with DOS p(e) - |p| r ,0 < r < 1) has been studied by 
a number of authors p5Q . They found a critical value 
for the Kondo coupling constant below which the local 
moment decouples. This feature has been reconfirmed 
in our calculations. However, beyond the critical value 
they found the same Fermi liquid strong coupling fixed 
point as in the standard Kondo problem. To the con- 
trary, some dramatically different results are obtained in 
our studies. We find a MFL behavior in the mixed va- 
lence and nearly empty orbital regimes. Namely, the real 
part of the self-energy goes like luIhuj, while the imagi- 
nary part behaves like \u>\ as u> — > 0. We believe a new 
universality class for the strong coupling fixed point has 
been found. The discrepancy with earlier treatments is 
due to the fact that different limits are considered. In 
their case the occupation of the impurity is always one 
and there is a true localized energy level well below the 
Fermi energy. We, in contrast, are considering the oppo- 
site limit, when the hybridization is assumed to be large 
and the impurity energy levels can merge with the con- 
duction electrons. From the theoretical point of view, the 
appearance of MFL behavior in our model is fully under- 
standable. It is well-known that near the nodes of a dSC, 
a Dirac type spectrum appears and the standard dimen- 
sional analysis of the quantum field theory can be ap- 
plied @. The scaling dimensions of the Nambu spinors 
^(r, t) and <p(t) turn out to be —1 and in length units, 
respectively. Thus short-range interactions between the 
conduction electrons are irrelevant, while the hybridiza- 
tion term of the Anderson Hamiltonian is marginal and 
is responsible for this MFL behavior. It seems to us that 
the Dirac structure of the energy dispersion itself is the 
main reason behind the MFL for the strong coupling fixed 
point. 

Focus now on N (r, oj) in the spatial range < r < £. 
Here £ = Kvf/Aq is the coherence length of dSC state, 
and also the natural length unit of our model, while in 
high T c dSC state, £ is about 10A, or roughly 3 lat- 
tice spacings. In Fig. 3, the local DOS vs frequency is 
shown for r=0.07£ from the magnetic impurity along the 
directions of the gap maxima and the gap nodes. In ad- 
dition to usual V-shape structure, there are quasiparticle 
resonances near the Fermi energy, and the positions of 
these resonant peaks coincide with those of the impurity 
resonances u> — ±fl'. Along the directions of the gap 
maxima, there are two resonances below and above the 
Fermi energy, which are slightly asymmetric in the line 



shape. On the other hand, along the directions of the gap 
nodes, there is only one sharp resonance and the local 
DOS is entirely hole- like. As the impurity energy level e<i 
increases, the quasiparticle resonances become broader, 
exhibiting a similar dependence as the local resonance of 
the impurity || . 

We also calculate the spatial variation of the DOS of 
the conduction electrons. The DOS around the magnetic 
impurity at the resonance energies is displayed in Fig. 4a 
for +0' and in Fig. 4b for —CI' as a function of spatial 
variables for e^/Ao = —0.2 in a logarithmic intensity 
scale. The quasiparticle resonances induced by the mag- 
netic impurity are highly localized around the impurity, 
and the spatial oscillation of these resonant states is visi- 
ble. The largest amplitude of the quasiparticle resonance 
occurs at the neighborhood of the impurity, and the lo- 
cal electronic structures distinctly differ in Fig. 4a and 
Fig.4b. For ui = +Cl', the local DOS exhibits a four- 
fold symmetry along the directions of the gap nodes for 
all distances, consistent with the dSC of the conduction 
electrons. For to = —CI' , the local DOS is strongly en- 
hanced in the gap maxima directions at distances r -C £. 
Further away from the impurity it is confined to 

the neighborhood of the diagonal directions, leading to 
an eight-fold symmetry. 

The logarithmic correction to the real part of the 
self-energy is a very subtle effect to detect experimen- 
tally. However, its imaginary part, the inverse quasi- 
particle lifetime r _1 oc — (V /ed) 2 Im oc \u>\ can 
be checked by experiments directly. Amazingly, such 
a linear frequency dependence of the inverse quasipar- 
ticle lifetime has been observed in the recent angle- 
resolved photoemission experiments in optimally doped 
Bi 2 Sr2CaCu20s+5 along the nodal directions ]T^Jl5|| . Al- 
though these high T c cuprates are believed to contain a 
small number of intrinsic defects or implications and " im- 
purity scatterings" may lead to localization of quasipar- 
ticle states, it is not clear whether the Anderson impurity 
model embedded in dSC state is applicable in this case. 

To conclude, we have investigated the quantum mag- 
netic impurity effects in high-T c superconductors based 
on the Anderson model. We have found a new univer- 
sality class for the strong coupling fixed point for this 
type of models. We have made explicit predictions on 
the resonance states around the magnetic impurities to 
be compared with experiments. 

After submitting the manuscript, we received a preprint 
on optimally doped Bi2Sr2Ca(Cui_ 2; Ni a ;)08+ ( 5 from Dr. 
S. H. Pan Jl6[ ] in which a localized resonance above the 
Fermi energy has been reported in the DOS of the Ni 
impurity, and the spatial dependence of the conduction 
electron DOS at the resonant energies is in a reasonable 
agreement with our calculations. 
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Figures Captions 

Fig. 1. The ground-state phase diagram of the model. 
Three areas denoted by I, II, and III correspond to decou- 
pled local magnetic moment, mixed valence, and strong 
coupling regimes, respectively. 

Fig. 2. The DOS Ni mp (lj) of the magnetic impurity, 
(a) for lj < and (b) for lj > 0, with e d /A = -0.2, 
0.0, 0.2, and 0.4, denoted by solid, dashed, dotted, and 
dash-dotted lines, respectively. 

Fig. 3. The local DOS N (r, lj) of the conduction elec- 
trons for ed/ A = —0.2,0.0, and 0.2 (from top to bottom) 
in units of Np. Here r=0.07£ correspond to the largest 
amplitude of the quasiparticle resonance at the neighbor- 
hood of the impurity, (a) along the directions of the gap 
maxima and (b) along the directions of the gap nodes. 

Fig. 4. The spatial distributions of the conduction 
electron DOS around the impurity at the resonant ener- 
gies, (a) lj = CI' and (b) lj = —CI'. Here Cd/Ao = —0.2 
and a logarithmic intensity scale is used. The coherent 
length ( is about 10A, or roughly 3 lattice spacings in 
high T c cuprates. 



4 



